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We construct a large class of completely positive and trace preserving non-Markovian dynamical 
maps for an open quantum system. These maps arise from a piecewise dynamics characterized by 
a continuous time evolution interrupted by jumps, randomly distributed in time and described by a 
quantum channel. The state of the open system is shown to obey a closed evolution equation, given 
by a master equation with a memory kernel and a inhomogeneous term. The non-Markovianity of 
the obtained dynamics is explicitly assessed studying the behavior of the distinguishability of two 
different initial system's states with elapsing time. 
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Open quantum systems naturally arise in quantum me- 
chanics due to lack of isolation, and one of the basic 
difficulties in the field is the derivation of closed irre- 
versible evolution equations for the system only, taking 
into account the interaction with the environment ^[3]. 
fn particular an open issue is the characterization and 
study of memory effects described by these irreversible 
dynamics. An important class of dynamical evolutions 
is given by quantum dynamical semigroups, which by 
construction ensure complete positivity (CP) and have 
a number of attracting physical and mathematical fea- 
tures. The semigroup property ensures the existence of 
a closed evolution equation, known as master equation, 
whose general expression has been determined in the 70 's 
just thanks to the requirement of CP [H[S]. The opera- 
tors appearing in the the master equation can be easily 
linked to the microscopic events which characterize the 
dynamics. Moreover the exact solution can be expressed 
in terms of a Dyson expansion, which allows for a natural 
reading in terms of a piecewise dynamics consisting of a 
relaxing evolution interrupted by jumps. 

In this Letter we show how a similar construction can 
be exploited to obtain a large class of non-Markovian 
completely positive trace preserving (CPT) maps, still 
admitting closed evolution equations. The building 
blocks of this construction are a collection of time de- 
pendent maps, together with a waiting time distribution 
describing the random occurrence in time of interaction 
events described by a quantum channel. The operational 
construction provides a direct physical reading of the dif- 
ferent contributions to the dynamics. The resulting mas- 
ter equations exhibit an integral kernel which warrants 
CP of the solution, one of the crucial difficulties in look- 
ing for extensions of the Lindblad result [BHU] . 

Master equations. For a semigroup we have p(t) = 
$(i)iO, where the time evolution operator obeys the mas- 
ter equation d$(t)/dt = £&(t) and satisfies 

*(*i + t2) = *(*a)*(*i), Vti,t 2 >0. 

Introducing a self-adjoint operator H and operators L k 
that can be associated to microscopic interaction events, 



e.g. the exchange of an excitation between system and 
bath, the operator L called the generator takes the 
form [1 [S] Cp = Rp + pR) + Jp, where R = -iH - 
(1/2) J^k L\Lk, and the CP superoperator J reads 

Jp = Y J L kP L\. 

k 

Introducing further the superoperator lZ(t), which gives 
the semigroup obtained exponentiating the operator R 

TZ(t)p = e tR pe m \ 
the exact evolution can be written as the Dyson series 

1>(t)p = K{t)p+J2 dt n ... dti (1) 
xTZ(t - t n )JTZ(t n - t n -i) . . . jn{h) P . 

Here p denotes the reduced system state taken as ini- 
tial condition, and the result follows from the Schwinger 
formula |I0j granting in particular trace preservation. 
This solution can be naturally described as a sequence 
of jumps, corresponding to transformations induced by 
the CP map J ', distributed over an underlying relaxing 
evolution given by the semigroup lZ(t). This kind of dy- 
namics is universally accepted as Markovian. Indeed the 
fact that the state of the system at a time t\ + t^ only 
depends on its state at a previous time t\ expresses a 
feature that is naturally associated to lack of memory 
and therefore to Markovianity (M). In this sense also a 
collection of two time evolution maps $(£ + t, t) obeying 
the composition law 

$(ti + t 2 ,0) = 3>(t 1 + t 2 ,t 1 )$(tx,0), Vti,i a > 

where each map is CPT, embodies the same idea of inde- 
pendence from the states at previous times, and is there- 
fore taken as a natural criterion to assess or define M, 
known as divisibility [IIJ [12] . Most recently a novel idea 
has been put forward to characterize M, neither basing 
on a representation of the dynamics, nor on the notion 
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of memory as dependence on the previous states of the 
system, but rather on the notion of distinguishability of 
system's states, and on its behavior in the course of the 
dynamics, which calls for an involvement of the environ- 
ment and of correlations [T31 [T5] . It turns out that this 
criterion is satisfied by a dynamics characterized by di- 
visibility, but is in general less restrictive [T5Vll8j . 

Derivation from piecewise dynamics. We now build 
on these known results to construct a much wider class of 
time evolutions, which admit a natural reading in terms 
of a piecewise dynamics, with microscopic interaction 
events embedded in a continuous time dynamics. These 
dynamics obey closed evolution equations expressed by 
means of a master equation, possibly admitting an inho- 
mogeneous contribution, which keeps track of the initial 
condition. As a starting point we consider Eq. ([I]), replac- 
ing the semigroup 1Z(t) with a collection of time depen- 
dent CPT maps F{t), which describe the time evolution 
between jumps. The events taking place over the back- 
ground of the continuous time evolution are described by 
a CPT map £ , namely a quantum channel, and their dis- 
tribution in time is characterized by an arbitrary waiting 
time distribution, so that the number of events in time 
realizes a renewal process. In terms of these basic build- 
ing blocks one has, given an initial state p, a time evolved 
state given by 



A(t)p = p (t)F(t)p- 

X p n ini ■ 



oc 

E 

n=l 



At, 



dii 



,t 1 )F(t-t n )£...£T{t 1 )p. (2) 



Here p n (t;t n , . . . t ti) denotes the exclusive probability 
density for the realization of n events up to time t, at 
given times t\, . . . ,t n , with no events in between. This 
probability density for a renewal process reads 

p n (t; t n , ...,h) = f(t- t n ) . . . f(t 2 - ti) S (ti), (3) 

with f(t) a waiting time distribution, i.e. a distribution 
function over the positive reals, and g(t) = 1 — J Q dr/(r) 
its associated survival probability, expressing the proba- 
bility that no jump has taken place up to time t. Thanks 
to CPT of the maps £ and JF(t) the obtained dynamics is 
indeed well defined. CP is warranted by stability of the 
positive cone of CP maps under composition. Regarding 
trace preservation, due to Eq. ^ for a renewal process 
the probability p k (t) to have k counts up to time t obeys 



p k (t)= / Arf{t-r)p k ^{r), (4) 
Jo 

with po(t) = g(t). Iterating this identity one obtains 
TrA{t)p — Y^k=oPk{t)P = P- The constructed collection 
of CPT time evolutions A(t) are functionals of T(t), f(t) 
and £, and allows for a simple operational interpreta- 
tion in terms of the random action of a fixed quantum 
channel over a given dynamics, not necessarily obeying a 
semigroup composition law. 



Laplace transform and master equation. We now ob- 
serve that, according to its definition Eq. ([2]), the map 
A(i) obeys the integral equation 

A(i) = g(t)T(t) + f dr/(t - r)F{t - t)£A(t), (5) 
Jo 

which in Laplace transform, here denoted by a hat, sim- 
ply reads 



A(u) = gF{u) + fT(u)£A(u). 



(6) 



Starting from this expression, as described in the Supple- 
mental Material [H] one finally obtains the closed master 
equation 



P (t) = / drK{t - r)£p(r) + J(t)p(0), (7) 



dr 



with kernel and inhomogeneous term given by 



IC(t) = -[f(t)T(t)}+f(0)d(t) l{t) = -[g{t)T(t)\. (8) 



dt 



dt l 



This is the main result of our Letter. We stress the fact 
that the map A(i), solution of Eq. (j5j), or equivalently 
Eq. |7b, is CPT by construction. It can be obtained as 
the inverse Laplace transform of the solution of Eq. ^ 



A(u) = 1 - fF{u)£ gF{u). 



(9) 



This identity provides a compact general expression of 
the Laplace transform of the exact solution, in terms of 
the transform of the elementary maps determining the 
time evolution. Note that the result has been obtained 
without making any restrictive assumption on the dimen- 
sionality of the Hilbert space of the system. 

Limiting expressions Before considering the non- 
Markovianity (NM) of the class of master equations in- 
troduced above in view of the recently proposed criteria 
[TT1 [T5] , we to point to some special cases already con- 
sidered in the literature. Firstly a quantum dynamical 
semigroup is recovered if F{t) — > e tc , with C in Lind- 
blad form, and £ — > 1, independently of the waiting time 
distribution f{t). Indeed, the solution given by Eq. Q 
thanks to the properties of the Laplace transform with re- 
spect to shifts now reads A{u) = X^feLo 5( u ™£)/ fc ( M ~ £)> 
and therefore, also using Pk(u) = g(u)f k (u), which fol- 
lows from Eq. Q, we have p(t) — e tc p. More generally, 
for a non trivial CPT map £ rearranging terms one ob- 
tains [TS] 



p(t)=Cp(t)+ / dTk{t-T)c^- T ^ c [£ -1] P (t), (10) 



dt' 



where the (D-number kernel reads k{u) — f{u)/g{u). 
This equation has been previously considered for the 
special case of a Lindblad generator given by a simple 
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commutator, pointing to a possible microscopic deriva- 
tion [3 [20]. For a vanishing Lindblad generator one 
has in particular p(t) = J2kLoPk{t)£ k P> a class of non- 
Markovian evolutions studied in [TH1 [H] . 

If we allow for a generic CPT map J 7 it), but do con- 
sider the events as a reset of the continuous time dynam- 
ics described by Tit), so that £ — > 1, we end up with 



-^)= J &Tf{t-T)F{t-T)p(T)+g{t)F{t)p, (11) 

which for the case of a memoryless waiting time of expo- 
nential type, fit) — Te~ rt , so that git) = e _r ', recovers 
the result recently obtained relying on a collisional model 
assuming collisions with independent ancillas |22j . 

Non-Markovianity. We now study the NM of the dy- 
namics described by the master equation Eq. Indeed, 
despite the fact that the considered master equation can 
include more general situations than a semigroup dynam- 
ics generated by a Lindblad operator, the degree of NM of 
the obtained dynamics is still to be ascertained. To this 
aim we will make reference to the definition of NM associ- 
ated to the idea of revival of distinguishability among dif- 
ferent states advocated in [T31 Q3] , considering the trace 
distance as a natural quantifier of distinguishability. As 
it has been shown, this criterion is more stringent than 
the violation of divisibility in terms of CP maps [T5HT8] . 
As a result, if we detect NM by using the notion of dis- 
tinguishability, we know that the considered dynamics 
is non-Markovian also from the divisibility point of view. 
We recall that the trace distance between two states p\ it) 
and P2(t) is given by the trace norm of their difference 
D(pi(i),j0 2 (i)) = 5 II Pi C*) — P2 (*) that is the sum of the 
modulus of the eigenvalues of their difference. It takes 
values between zero and one and can be interpreted as a 
measure of the distinguishability among states. In partic- 
ular, relying on the fact that the trace distance is a con- 
traction with respect to the action of a CPT map, M of 
the map is identified with the monotonic decrease in time 
of the trace distance among any couple of possible initial 
states. NM is then detected whenever the time deriva- 
tive of the trace distance grows at a certain time t, for 
at least a couple of initial states, i.e. D{p\{t), p%it)) > 0. 
In order to highlight this behavior, let us make specific 
choices for the system and the different maps and func- 
tions determining the time evolution A(t). We therefore 
consider the Hilbert space (D 2 , and take as CPT map £ a 
Pauli channel £ip — oipoi, with i = 0, x, y, z and do = 1 . 
We further take as waiting time distribution /(t) a con- 
volution of exponential distributions. These waiting time 
distributions bring with themselves a natural time scale 
given by the mean waiting time. Finally we have the 
freedom to consider a collection of time dependent CPT 
maps. The latter also have an intrinsic time scale, and 
the interplay between the two time scales plays an im- 
portant role in the characterization of NM. To this aim 
we will analyze two situations, corresponding to different 



physical implementations. As a first example we take 
a map !Fd{t) only affecting coherences, which according 
to the trace distance criterion by itself always describes a 
non-Markovian dynamics. As a complementary situation 
we will deal with a time evolution J- + it) which itself ad- 
mits both a Markovian and a non-Markovian limit, and 
affects all components of the statistical operator. 

Examples. We first consider a dephasing dynamics 
J^dit), which multiplies the off-diagonal matrix elements 
of the statistical operator by the function Working 
in (D 2 it is convenient to represent statistical operators 
through their coefficients on the linear basis {ctj/v2}, so 
that maps can be represented as matrices [531 122 • This 
dephasing map in particular is represented as a diago- 
nal matrix Fd(t) — diag(l,D(t),Dit),l), and the same 
holds for the Pauli maps which take the general form 
E = diag(l,e x ,s y ,e z ), with Si — ±1, the sign depending 
on the specific choice of map. Relying on Eq. ([£]), these 
expressions after some algebra [T5] lead to the following 
compact result for the time evolution map 



Adit) = diag(l,X(t),Y(t),Z(t)). 



(12) 



For the expression of the time dependent functions ap- 
pearing in the evolution map we consider the functional 



Lf [M] (u) 



gMju) 

l±fMiuY 



(13) 



where M denotes an arbitrary function of time. X(t) 
and F(t) are then given by one of the functions d±(t) = 
[D] it), depending on the value of the while Z{t) 
is given by either the identity or the function g(t) = 
Yl^=oP2n(t) - X^=oP2ri+i(£), which gives the difference 
between the probability to have an even and an odd num- 
ber of jumps. Given the explicit expression of the map, 
one can calculate the time derivative of the trace distance 
among two different initial states, which shows in partic- 
ular that one has NM whenever the modulus of one of the 
functions d±it) or qit) grows, as discussed in the Supple- 
mental Material [15]. This case is depicted in Fig. [lja), 
considering a dephasing map Z?(i) = cos(Ai). In this 
case the dynamics given by Tdit) alone never allows for 
a Markovian description. Here the rate A sets the nat- 
ural time scale for this contribution to the dynamics, to 
be compared with the time scale 1/r given by the mean 
waiting time associated to the waiting time distribution 
fit). As it appears in Fig. [lja), if T/A 3> 1, so that sub- 
sequent events are very close in time, the contribution to 
NM due to Tdit) is suppressed, since on a short enough 
time any time evolution map is Markovian. 

As a further example we consider the dynamical map 
J-+it), affecting both populations and coherences, that 
arises considering the interaction of a two-level system 
with a bosonic field in the vacuum state pQ. The map 
is characterized by the function Git), depending on the 



4 




FIG. 1. (Color online) (a) Modulus of the functions d_(t) and 
q(t) for a dephasing dynamics described by D(t) = cos(Af), 
and waiting time given by the convolution of three equal ex- 
ponentials re _r *. The growth of any of these quantities, as 
discussed in the Supplemental Material [T5j, provides a direct 
signature of NM of the time evolution map Ad(t). The quan- 
tities are plotted as a function of Xt and F/X. The semitrans- 
parent surface corresponds to d-(t), while the meshed surface 
represents q(t). It immediately appears that for growing ratio 
F/X, determining the relation between the time scales inher- 
ent in J~d{t) and f(t), the oscillations in d-(t) are suppressed. 
The NM is then only detected by q(t), arising due to the 
action of the map £ x which describes the events, in this case 
spin flips, in between the continuous time evolution J-d(t). (b) 
Modulus of g-(t) and h+{t), here for a continuous dynamics 
J-+(t) involving both populations and coherences, and waiting 
time corresponding to the convolution of two equal exponen- 
tial distributions. The ratio 7/ A appearing in the function 
G(t) given by Eq. ( 16 1 is set equal to 3, corresponding to NM 
of J~+(t) alone. Again the growth of the modulus of any of 
these functions warrants NM. The function g~(t) given by the 
semitransparent surface only detects NM for small F, while 
h+ (t) shows that NM also takes place for frequent events, that 
is large F/X, even though it is confined to shorter and shorter 
times. 



spectral density of the environment, and in matrix form 
reads 

F + (t)=diag(l,G(t),G(t),\G(t)\ 2 ) + B(\G(t)\ 2 -l), (14) 

where B(x) denotes the 4x4 matrix with entry x in the 
bottom left corner as the only non zero element. Ex- 
ploiting Eq. ([9]) we can obtain the expression of the time 
evolution map 19] 

A+(t) = diag(l, X(t),Y(t), Z(t)) + B{W{t)), (15) 

where now X(t) and Y(t) take the expressions 
g±(t) = [G] (t), while Z(t) corresponds to h±(t) = 
L f[\G\ 2 ]{t). The function W(t), determined by /(*) and 
|G(i)| 2 , does not affect the trace distance, since it corre- 
sponds to a fixed translation of the state [25 . A typical 
expression of G(t) is given by 

G(t) = e~ xt/2 [cosh( 7 i/2) + (A/7) srah(7*/2)] , (16) 

where 7 = \/X 2 — 27A, and has the interesting feature 
that for 7/A < 1/2 the map J-+{t) itself is Markovian, 



while for 7/A above this threshold one has NM [IS]. The 
NM of the ensuing overall dynamics A(t) is considered 
in Fig. |ljb) , where we have plotted the modulus of the 
functions g-(t) and h + (t) for G{t) as in Eq. (16). Again 
the growth of the modulus of any of these functions is a 
witness of NM. It appears indeed that for a wide range of 
parameters the dynamics is non-Markovian, yet the NM 
is actually the result of an interplay of the features of 
all the three elements determining the dynamics, namely 
T + {t), £ and f{t). Indeed for values of the ratio 7/A 
such that F+(t) itself is non-Markovian, the dynamics 
A(i) might still be Markovian, if the ratio T/X of the 
time scales associated to J-+(t) and f(t) is high enough. 
On the contrary, even a Markovian J- + (t) can give rise to 
a non-Markovian dynamics because of the action of the 
map £ in between the continuous time evolutions, and of 
the distribution in time of these events. 

Conclusions. We have obtained a large set of closed 
non-Markovian master equations starting from a piece- 
wise dynamics described by a continuous time evolution 
interrupted by random jumps. The solution of these 
equations is warranted to be a CPT map. These master 
equations involve both a memory kernel and a inhomo- 
geneous term. The basic ingredients in the construction 
are a collection of time dependent maps, together with 
a waiting time distribution describing the random occur- 
rence of events characterized by a quantum channel. We 
have considered the connection of this result with the 
standard expression of quantum dynamical semigroups, 
as well as more recent examples of non-Markovian mas- 
ter equations obtained starting from microscopic models. 
In particular, we have certified the NM of the obtained 
time evolution by studying the behavior in time of the 
distinguishability between two different initial states, as 
quantified by the trace distance. Finally, the operational 
interpretation of the structure of these master equations 
paves the way for their use in concrete applications. 
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SUPPLEMENTAL MATERIAL 

In this Supplemental material we provide technical de- 
tails on the derivation of equations and properties dis- 
cussed in the main text of the paper. 



Derivation of Eq. 

We here derive the closed master equation obeyed by 
the statistical operator pit). Given that A(i) obeys the 
integral equation Eq. ([5]), as considered in the main text 
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in Laplace transform the equation for A(u) reads 
A(u) = $F{u) + fT{u)£k{u), 

with the Laplace transform defined as usual and denoted 
by a hat h(u) — dte~ ut h(t), so that multiplying by 
u and subtracting the identity operator from both sides, 
at the same adding and subtracting the term f(0)£k(u) 
at the l.h.s. one comes to 



uk{u) 



ugF(u) — 1 
+/(0)£A(u), 



ufT(u)-f(0) £A(u 



so that recalling that the Laplace transform of the deriva- 
tive of a function h(t) is given by uh(u) — 1, and using 
J-{0) = 1, one obtains 



dt 



A(i) = / dr 
Jo 



d(t-r 



-f{t-T)F{t-T)£k{r) 



d 



-f(0)£ A(t) + - fo(t).F(t)] 



According to the relation = A(t)p(0) and using the 
identifications Eq. ^ one finally comes to the master 
equation Eq. 



Derivation of Eq. and Eq. J77 



In order to derive the master equation Eq. ( 10 1 we start 
from Eq. (JgJ) and take J-{t) to be of exponential form e tc 
with C a Lindblad generator. Thanks to the behavior of 
the Laplace transform with respect to shifts one thus has 
gf(u) = g(u — £), and similarly for f^Fiu), so that for 
the Laplace transform of the statistical operator p(u) we 
obtain 

p{u) = g(u - £)p(0) + f(u - C)£p(u). (17) 

To proceed further we note that from the relation be- 
tween waiting time distribution and survival probability 



g(t) = 1 - / dr/(r), 
Jo 



(18) 



one has 



1 - /(«) 



and therefore introducing the function 



k{u) = 



(19) 



also 



1 



k(u - C) 



C. 



g(u - C) 

We note that the function k(t) naturally appears as 
memory kernel in the description of continuos time ran- 
dom walks [53]. Dividing Eq. p7]) by g(u — C) and 



using Eq. ( |19[ ) one thus obtains, subtracting a term 
k(u — C)p(u) from both sides 

up(u) - p{0) = Cp(u) + k{u -£)[£- 1] p(u), (20) 

and finally taking the inverse Laplace transform, exploit- 
ing again the property of the Laplace transform with re- 
spect to shifts, the master equation 



dt' 



p{t)=Cp{t)+ / &rk{t - t)c^-^ c [£ -1]p(t). 



For the case of Eq. (11 1 we start from Eq. ^ and again 
multiply by u and subtract the identity from both sides, 
so that suitably rearranging terms and taking £ — > 1 wc 
have 



uA(u) - 1 = ugT(u) - 1 + fF(u) uk(u) - 1 
+/>(«), 

leading to 
d 



dt A(t) = ^ dr/(t-r)^(t-T)A(T)+/(t)J-(t) 



di 



[ff(OT<)] 



and further exploiting the relation g(t) = —f(t) follow- 



11) 



ing from Eq. ( 18 ) one finally obtains the desired master 
equation Eq. 



A(t)= / drf(t-T)T(t-r)p(r)+g(t)T(t). 



dt 



Derivation of the map Ad(t) 

We now derive the time evolution map Ad(t) for a 
dephasing dynamics, which only affects the off-diagonal 
matrix elements of the statistical operator of the sys- 
tem, multiplying them by a function D(t), taken in 
the example to be cos(Xt). Any statistical operator on 
(D 2 can be represented by a vector in the linear ba- 
sis ^75 a x, ~^ a yi 01 'thonormal according to 
the Hilbert-Schmidt scalar product, so that maps can be 
identified with suitable 4x4 matrices O [24]. The 
dephasing map Fd{t) in this basis acts as the diagonal 
matrix 

F d (t)=diag(l,D{t),D(t),l), 

while the Pauli map £ x corresponds to the diagonal ma- 
trix 

E x = diag(l, 1, -1, -1). 

Starting from this result we have that the Laplace 
transform of the operator /(^^(t) can be written 
as diag(f(u),fD(u),fD(u),f(u)), and similarly for 



() 



g(t)Fd{t). Thanks to the closure of the algebra of diag- 
onal matrices Ad(u) itself turns out to be diagonal, and 
according to Eq. ([9| reads 



its derivative can be explicitly calculated. Using as in 



Arf(u) = diag | 
Upon defining 



f fD(u) fD(u) 11 -f(u) 



u l-fD(u) l + fD(u) ul + /(u) < 
gT){u) 



d ± {t) = Lj [D] (t) = 



1 ± fD(u) 



as in Eq. (131, as well as 



1 1 - f(u) 

u 1 + f(u) ' 



which according to the relation pk (u) — g(u)f k (u), which 
follows from Eq. j4]), is the Laplace transform of the 
quantity q(t) = £ n=0 P2n(i) - J2^=oP^+i(t), we finally 
obtain 

A rf (i) = diag{l,d-(t),d+(t),q(t)), 



which provides the explicit expression of Eq. ( 12 ) when 



the Pauli channel is given by £ x . Similar results apply for 
the other Pauli channels. The modulus of the functions 
d-{t) and q(t) is plotted in Fig. [IJa), since it provides 
evidence for NM of the dynamics, as discussed in the 
next paragraph. 



Non-Markovianity of the time evolution map 

We here apply the trace distance criterion for the de- 
tection of NM to the dynamics described by the map 
Arf(t), and similar conclusions hold for A + (i). As dis- 
cussed in the main text, according to this criterion NM 
is associated to the growth of the distinguishability in 
time, as quantified by the trace distance, of two distinct 
initial states. Given two initial states pi{0) and p 2 (0) one 
monitors their trace distance in time, as given by 

D( Pl (t), P2 (t)) = \\\pi{t) - P 2(i)\\i 

= \\\A d (t)(p 1 (0)-p 2 (0))\\ 1 , 

and the map is said to be non-Markovian if there exist 
a couple of initial states and a point in time such that 
their distinguishability grows, i.e. 



d 

dl 



D(p 1 (t), P2 (t))>0. 



For the case at hand, setting A p for the difference in the 
populations of the two initial statistical operators, as well 
as A c for the difference in the coherences, that is the off- 
diagonal matrix element, for a map diagonal in the basis 
used to represent states as vectors the trace distance and 



Eq. ( 12 ) the notation 

A d (t)=diag(l,X(t),Y(t),Z(t)) 

we have 

D( Pl (t),p 2 (t)) = 



A 2 p Z 2 (t) + Re 2 A c X 2 (t) + Im 2 A c Y 2 {t), 



and therefore 
d 



D( Pl (t),p 2 (t)) = 
1 A 2 p ±Z 2 (t) + fc 2 A4l 2 (t) + Im 2 A c Ay2( f ) 
A2Z 2 (i) + Re 2 A c X 2 {t) + Im 2 A c Y 2 (t) 



dt 



so that one has growth of the trace distance if the mod- 
ulus of any of the functions X(t), Y(t) or Z(t) grows. 



Derivation of the map A+ (t) 

We now consider as Pauli map £ z , and introduce a 
continuous time dynamics determined by a map F+{t) 
which in the above introduced basis for the operators in 
(D 2 is expressed as in Eq. (14) by the matrix 



F+{t) = diag(l, G(t), G(t), \G(t)\ 2 ) + B(\G(t)\ 2 - 1), 

where as discussed in the main text the matrix B(x) has 
the only non zero entry x in the bottom left corner. The 
calculations closely follow those performed for Ad(t). In 
particular thanks to the closure of the algebra of matrices 
with non zero entries only on the diagonal and in the 
bottom left corner, which are such that the inverse if it 
exists still is in the algebra, relying on Eq. ^ we obtain 



A + (u) = diag 



1 fG(u) fG(u) f\G\ 2 (u) 



1 + l + fG(u) i_/| G | 2 ( u)j 

+B(W(u)), 

where 

2/(«) - 1 + u (7\G?(u) - f\G?(u) 

W(u) = V 

l + f\G\ 2 (u) 

According to the definitions given in the main text below 



Eq. ( 15 ) we finally arrive to 



A+(<) = diag(l, h+(t), h+(t),g-(t)) + B(W(t)). 
Considering a map of the form 

A+(t) = diag(l, X(t),Y(t), Z(t)) + B{W(t)) 



7 



one immediately sees that the term W(t) provides a con- 
tribution to the matrix elements of the statistical oper- 
ator which is independent of the initial state, so that it 
does not affect the behavior of the trace distance. As a re- 
sult also in this case the M or NM of the map A + (t) does 
depend on the behavior of the modulus of the time depen- 
dent functions appearing on the diagonal of the matrix 
representation of A+(i). The latter is plotted in Fig.[]Jb). 
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